Oxygen Phonon Branches in Detwinned 
YBa2Cu307 

D. Reznik^'^, L. Pintschovius^ , W. Reichardt^, Y. Endoh'^, 
H. Hiraka^, J. M. Tranquada"^, S. Tajima^, H. Uchiyama^, and T. Masui^ 

^ Forschungszentrum Karlsruhe, Institut fur Festkorperphysik, 
Postfach 3640, D-76021 Karlsruhe, Germany 
^Laboratoire Leon Brillouin, C.E.A./C.N.R.S., F-91191-Gif-sur-YveUe 

CEDEX, France 

^Institute for Material Research, Tohoku University, Katahira, Aoba-ku, 
Sendai, 9808577, Japan 
'^Physics Department, Brookhaven National Laboratory, 
Upton, NY 11973 

^Superconductivity Research Laboratory, ISTEC, Shinonome, Koutu-ku, 
Tokyo, 135-0062, Japan 



We report results of inelastic neutron scattering measurements of phonon 
dispersions on a detwinned sample of YBa2Cu30^ and compare them with 
model calculations. Plane oxygen bond stretching phonon branches disperse 
steeply downwards from the zone center in both the a and the b direction in- 
dicating a strong electron-phonon coupling. Half way to the zone boundary, 
the phonon peaks become ill-defined but we see no need to invoke unit cell 
doubling or charge stripe formation: lattice dynamical shell model calcula- 
tions predict such behavior as a result of branch anticros sings. There were 
no observable superconductivity-related temperature effects on selected plane 
oxygen bond stretching modes measured on a twinned sample. 

PACS numbers: 78.70.Nx , 74.72.Bk, 74.25. Kc. 



Charge degrees of freedom associated with lattice vibrations and/or 
distortions may be important to our understanding of many physical prop- 
erties of cuprate superconductors. Longitudinal vibrations of planar oxy- 
gens soften strongly with doping away from the zone center in all cuprates 
investigated so faroa, which has been interpreted as a signature of strong 
electron-phonon coupling. Up to now twinning has made it difficult to ob- 
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Fig. 1. Measured (a) and resolution-corrected calculated (b) neutron scat- 
tering intensity at Q = (4-X,l,0) r.l.u. Background, linear in energy and 
independent of Q, has been subtracted from the experimental data. 



tain a clear experimental picture of this behavior in the middle of the zone, 
where phonon dispersions are steep. To clear up twinning-related ambigui- 
ties, we measured phonon dispersions by inelastic neutron scattering on an 
untwined sample of YBa2Cu307. High-resolution measurements on a large 
twinned sample confirmed these results. 

The untwinned sample consisted of close to 30 detwinned single crystals 
with the mosaic spread of 3° and the total volume of 0.8cm^. The twinned 
sample consisted of 3 single crystals of combined volume of 1.5cm^ with the 
mosaic spread of 2.2°. The experiments were performed on the triple-axis 
spectrometer IT at the Orphee reactor using doubly focusing monochroma- 
tor (Culll and Cu220) and analyzer (PG002) crystals, a-axis and b-axis 
atomic vibrations of the Al symmetry were measured in the Brillouin zone 
adjacent to wavevectors, Q, of (4,1,0) and (-1,4,0) respectively. The data 
were fit with gaussian lineshapes on top of a background assumed to be 
linear in energy, w, and independent of Q. 

The starting point for our calculations was the common interaction po- 
tential modelo, which is quite successful in describing the phonon dispersions 
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Fig. 2. Measured (a) and resolution-corrected calculated (b) neutron scat- 
tering intensity at Q = (-1,4-X,0) r.l.u. Exactly the same background as in 
Fig. 1 has been subtracted from the experimental raw data. 



of a number of cuprates. YBa2Cu307 is treated as an ionic compound and 
the interatomic interactions are modelled as a sum of Coulomb forces and 
short range repulsive forces. The polarizability of the atoms is treated using 
the shell model formalism. For a metallic compound like YBa2Cu307 a term 
accounting for screening by free carriers is added. This model gives a decent 
description of the available data, which was further improved by tuning the 
parameters. Further tuning was done to account for the new data. 

Such a simple model cannot reproduce the anomalous softening of the 
bond-strctching phonons, and thus would not realistically simulate the ex- 
periment. In particular, it would miss the effect of the anomalous softening 
of the bond-stretching modes on other phonons: it leads to anti-crossings 
with branches of the same symmetry, which originally have a different po- 
larization, but hybridize with the bond-stretching branches if they are close 
in energy. Special terms were added in order to model such effects: 1) A 
negative breathing deformability lowers the energy of the planar breathing 
mode; 2) A term lowers the energies of the 'half-breathing' modes at the 
zone boundary of the (100) and - with a different parameter - of the (010) 
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direction; 3) Similar terms produce maximum softening half way to the zone 
boundary in the 100 and 010 directions, but have no effect both at the zone 
center and at the zone boundary. The physical meaning of the special terms 
(especially the last two) is not completely clear. We think that they mimic 
the effects of a strong electron-phonon coupling. The calculations also in- 
cluded a resolution correction in Q and to. 

The longitudinal plane oxygen vibrations with atomic displacements in 
the a and b directions have the zone center energies of 74 and 68 meV re- 
spectively (Figs. 1, 2). Their zone boundary energies are 58 and 55 meV 
respectively. As the a-axis polarized longitudinal phonon disperses down- 
ward from the zone center, the peak disappears almost entirely at qa=0.15 
and 0.2 and then reappears at higher q and lower energy (55meV). Along 
the b-axis, the peak similarly disappears and reappears somewhat closer to 
the zone center. The difference between the energies of the a and the b-axis 
vibrations is entirely accounted for by the difference in the bond lengths. 
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Fig. 3. Observed and calculated phonon frequencies for the Al branches. 
The lines and the dots represent calculated and experimentally measured 
phonon energies for q perpendicular (a) and parallel (b) to the chains. 
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These data seem to imply that the dispersions of the longitudinal oxygen 
vibrations are discontinuous and one may be tempted to conclude that the 
lattice dynamics in YBa2Cu307 are highly unconventional. Similar data 
have been explained by unit cell doubling or charge stripe formationill. 

Inspection of the figures shows that the final model reproduces both 
the experimental frequencies and the observed intensities quite well (Figs, 
lb and 2b). It predicts that an almost flat apical oxygen branch should 
cross the downward-dispersing longitudinal plane oxygen branch. The apical 
oxygen branch does not show up in figures 1 and 2 because of its small 
structure factor in the measured Brillouin zones (in fact, zero at T). Since 
the apical and plane oxygen branches have the same symmetry, they do 
not actually cross, but mix and repel each other lifting the degeneracy at 
the crossing point. This is the branch anticrossing phenomenon described 
above. In the process, the lower branch acquires plane oxygen character on 
the higher q side of the crossing point and becomes visible in the measured 
Brillouin zones. Simultaneously, the upper branch acquires mostly apical 
oxygen character and becomes much weaker. 

Figure 3 shows that there is a good agreement between the experimen- 
tally observed phonon peak positions and the calculated dispersion curves. 
Furthermore, the shell model predicts the structure factors that at least 
qualitatively account for experimentally observed phonon intensities (Figs. 
1, 2). Thus we see no need to invoke unconventional theories such as charge 
stripes or unit cell doublings to describe lattice dynamics of YBa2Cu307. 

Finally, we found no appreciable temperature dependence of longitudi- 
nal oxygen vibrations measured on a twinned sample other than slight soft- 
ening with increasing temperature consistent with increasing anharmonicity. 
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